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We investigate the space-time structure of the classical gluon field produced in an ul- 
trarelativistic collision between color charges. The classical solution which was computed 
^ ' previously in a perturbative approach is shown to become unstable on account of the non- 

, Abelian self-interaction neglected in the perturbative solution scheme. The time scale for 

' growth of the instabilities is found to be of the order of the distance between the colliding 

^5 [ color charges. We argue that these instabilities will eventually lead to thermalization of 

gluons produced in an ultrarelativistic collision between heavy nuclei. The rate of thermal- 
ization is estimated to be of order g^fi, where g is the strong coupling constant and fi^ the 
transverse color charge density of an ultrarelativistic nucleus. 
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^ ■ I. INTRODUCTION 

c ■ 

Ultrarelativistic heavy-ion collisions at RHIC and LHC aim towards an understanding of nuclear matter 
k>( ' properties under the extreme conditions of high energy density and small (net) baryon number density 

r> I [Q. Lattice simulations of quantum chromodynamics (QCD) predict the occurrence of a (phase) transition 

■ between ordinary hadronic matter and a phase of restored chiral symmetry, where quarks and gluons are 

deconfined, at a temperature Tc 150 MeV [^j. Such a phase, commonly termed quark-gluon plasma 
(QGP) 1^, is believed to have existed in the early universe. The relative smallness of Tc (in comparison 
to the electroweak transition temperature) will, for the first time, enable us to study a phase transition 
predicted by one of the fundamental theories of matter under laboratory conditions. 

In order to assess whether the QGP can be re-created in ultrarelativistic heavy-ion collisions, high-energy 
nuclear theorists are challenged to clarify whether highly excited nuclear matter can actually equilibrate 
on time scales that are short enough such that (local) temperatures still exceed Tc- If matter equilibrates 
at later times, the system will already have decayed into hadrons. Then, possible experimental signatures 
of the high-energy-density phase might be washed out by more complicated and thus less well understood 
non-equilibrium phenomena. 

The issue of equilibration requires understanding of the kinetic processes in the pre-equilibrium phase of 
an ultrarelativistic nuclear collision. These, in turn, can only be understood knowing the initial conditions 
for the collision in the first place. Recently, McLerran and Venugopalan Q have made considerable progress 
in that direction. In their model of an ultrarelativistic nucleus, the valence partons act as a classical color 
source for small-a; gluons with transverse momenta ^qqd <C fc^ ^ where Aqcd — 200 MeV is the QCD 
scale parameter and /i^ the transverse color charge density of an ultrarelativistic nucleus. The reason why 
the color source can be treated classically is the following: gluons with sufficiently small x do not resolve 
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the nuclear substructure, the (Lorentz-contracted) nucleus appears coherent in longitudinal direction. In 
transverse direction, color is of course confined on scales larger than the nuclcon's size. The transverse 
momenta of the respective gluons are, however, sufficiently high to coherently probe the colored interior of 
individual colorless nucleons. Therefore, all valence partons throughout the longitudinal extension of the 
nucleus and inside a transverse area ~ add up to the color source "seen" by these gluons. Obviously, 
the transverse color charge density /i^ grows like A^/^ and thus becomes large for heavy nuclei. The color 
charge in a given transverse area is therefore in a high-dimensional representation of the color algebra, and 
thus essentially classical. 

The original model of McLerran and Vcnugopalan has recently been extended to include contributions of 
the sea partons |^,^. Furthermore, the color charge density /i^ has been calculated in a more microscopic 
model of the nucleus and limitations of the classical approach have been pointed out Also, the classical 
gluon field produced in a collision of two such nuclei has been computed in a perturbative scheme ||9|,|l0|] and 
has been shown to agree with perturbative gluon bremsstrahlung as calculated in The classical field 



was shown to have a quantum interpretation as a coherent state 1 12 . The relationship between the classical 
result and mini-jet production in the conventional coUinear factorization approach JTst as well as from a 
BFKL-ladder has been discussed in Q. 

In this work, we aim to clarify the space-time structure of the classical gluon field computed in |lO[ and 
investigate its stability. It is well known that, in contrast to Abelian theories, non-Abelian gauge fields are in 
general classically unstable, due to their spin-field interaction with the corresponding gyromagnetic ratio for 
the gluon (see for instance A well-known case of this instability is the decay of a constant color-magnetic 

field The dynamical instabilities have been shown to generate classical chaos [Q. Furthermore, it was 
demonstrated [Q that they lead to inelastic final states in high-energy collisions of wave packets. Here we 
shall show that the gluon field produced in ultrarelativistic nuclear collisions has at least one unstable mode 
throughout the evolution of the system. The growth rate of the instability is determined to be of order g^fj,. 
The unstable mode originates from the self-interactions described by the non-Abelian terms in the Yang- 
Mills equations, which have been neglected in the perturbative solution HJl^. Since these self- interactions 
may, under favorable conditions, lead to thermalization, the rate of growth of the instability also sets the 
time scale for thermalization of the produced (classical) gluon field. 

The remainder of this work is organized as follows. In Section II we briefly summarize the derivation of the 
(perturbative) classical gluon field |l^] . In Section III we show the detailed space-time structure of this field 
for a collision between two individual color charges. The complete solution in the case of a nuclear collision 
is simply a linear superposition of such color fields. Section IV has three parts. In the first, we give a general 
outline of the linear stability analysis which will be employed to investigate the stability of the perturbative 
solution. In the second part, we apply this scheme to re-derive the decay of a constant color-magnetic field. 
Finally, in the third part, we investigate the stability of the classical gluon field discussed in Sections II and 
III. In Section V we summarize our results and draw conclusions for the thermalization rate. 

Our units are h — c = I, and the metric tensor is g^'^ — diag(-f, — , — , — ). Light-cone coordinates are 
defined as a± = (a°±a^)/V2, dzp = d/dx±. The notation for 3-vectors is a = (a^, a^, a^) and for transverse 
vectors a = (a^, o^). 



II. CLASSICAL GLUON RADIATION IN ULTRARELATIVISTIC NUCLEAR COLLISIONS 

In this section we briefiy outline the derivation of the classical gluon field produced in an ultrarelativistic 
nuclear collision (for details, see Ref. |l§). In covariant (Lorentz) gauge, d^A^^ — 0, the classical Yang-Mills 
equations, 

D^F'^'' - r , (1) 

where I?^ = — ig [ , • ] , F^'' = d^^A'^ — d'^ A^^ — ig [A^ , A'^], and J'^ is a classical source current, can 
be cast into the form 

nA" = + ig [ A^ , d'^A" +F''''] , (2) 

where □ is the d'Alembertian operator. In this form, the equations can be solved perturbatively order by 
order in the strong coupling constant, g. We expand 
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fe=0 fc=0 

where A^jf^ , J^"'' are the contributions of order to the gluon field and source current, respectively. (Eq. (^) 
takes into account that only odd integers of g occur in this expansion.) Then, the lowest and next-to- lowest 
order solutions obey 

□41) = J(i) ^ J(i) , (4a) 
= ) +^g[A^'>. a. A W + Fii) ] ^ ) . (4b) 



These equations are linear to each successive order in g and can therefore be solved with the method of 
Green's functions: 

42/c+i) ^ f ^4^/ ^^(^ _ j(2fc+i) ^ A: = 0, 1 . (5) 



The retarded Green's function reads in coordinate space [g^ : 

Gr{x)^^e{t)5{x') . (6) 

The explicit solution requires the specification of initial conditions. Here we consider a collision of two nuclei 
with mass numbers Ai, A2, moving towards each other with ultrarelativistic velocities, wi^2 — il, along the 
2;-axis. The nuclei are taken as ensembles of nucleons IfTlEOI, cf. Fig. 1. 




FIG. 1. The nuclear collision as envisaged here. 



In order to simplify the color algebra we assume that each "nucleon" consists of a quark-antiquark pair. 
These valence quarks and antiquarks are confined inside the nucleons (visualized as spheres of equal radius 
in the rest frame of each nucleus). Nucleons inside the nucleus and valence charges inside the nucleons are 
assumed to be "frozen", i.e., they have definite light-cone and transverse coordinates. We assume the charges 
to move along recoilless trajectories, therefore, their coordinates will not change throughout the calculation. 
We label the coordinates of the quarks in nucleus 1 by X-i, x^, i — I, . . . ,Ai, and those of nucleus 2 by 
y+j, y j = I, ■ ■ ■ , A2. Antiquark coordinates follow this notation with an additional prime. 

Then, the lowest-order classical current is a sum of the currents for each individual nucleus, as given in 
Eq. (3) of Re f. ||lOt| . The lowest order solution, A'j}\ and the associated field strength tensor, Pj^]) , were 
given in Ref. ]lO||7Eqs. (4,5). This solution is identical to the one of the corresponding Abelian problem, i.e.. 
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the nuclei pass through each other without interacting. The fields generated by the valence charges simply 
superpose and no gluons are radiated. 

(3) 

With this lowest-order solution and the classical current to next-to-lowest order in g, J/i , one can compute 
the next-to-lowest order solution, A^f^ . The classical current jj^'^ was obtained in [0 from covariant current 
conservation, D^J'^ = 0, and the assumption of recoilless trajectories, see Eq. (8) of The solution 

fr>\ 

; a;^, Ui) for a collision between quark i of nucleus 1 and quark j of nucleus 2 was found to be (cf. Eqs. 
(21-24) in 0]) 
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(7b) 



(7c) 



Here, A° = f^ix^lTl) (Tj), where (T^ ), (Tj) are color matrices which represent the color charge of the 
quarks in the color space of nucleon i of nucleus 1 and nucleon j of nucleus 2, while /"'"^ are the structure 
constants of SU(-/Vc). A is an infrared cut-off which acts as a gauge parameter, C>(<) = '[nax (min)(T, |^— xj), 
ri>(<) — max (min)(r, \x—y^\), where r = y^2(a;_ — X-i){y+ — y+j) is the proper time of the collision between 
quarks i and j. Finally, r = {(,<V<)/{(,>V>), cosa = {x-x^ ■ {x ~ y.) / {\x - x^\ \x~y.\), and = d/dx^. 
The complete solution for the nuclear collision is 

4"°(^) = E [4'^'^(^;^.,y,)-4'^n2:;x^,2/,)-4''°(^;^-y^-) + 4'^"(^;^^,?/;-)] ■ («) 



The antiquarks (with the primed coordinates) have the opposite color charge, —{Tj'), —{Tj), which explains 
the relative signs in Eq. (0) , and ensures color neutrality of each nucleon. 



III. SPACE-TIME STRUCTURE OF CLASSICAL GLUON RADIATION 

The classical gluon field (||) is a linear superposition of the fields generated in two-body collisions between 
the valence partons. In order to clarify its space-time structure it is therefore necessary to first consider 
the field produced in a single parton-parton collision, Eqs. ([7a|-|7^). For the sake of simplicity, we assume 
X-i = Xi = 0, and y+j = Q,y. = h= (6,0). Furthermore, we choose (by means of a suitable global color 

rotation) the orientation in color space such that A" = f"'^''{Tj'){T^) = AS""^. For times t > 0, the radiated 
gluon spectrum is boost-invariant along the 2:— axis |p|,p|pT|], therefore it is sufficient to consider the transverse 
plane at z = 0. It is furthermore convenient to choose the gauge parameter X = 1/b. Then, the non-zero 
components of the gauge field simplify to (we omit the superscript "(3)" and the arguments Xi, yj in the 
following) 
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A\{t,x,0) = -1- { 0(i - \x\) ln(r;>/6) + e{t - \x-b\) ln(C>/5) 

+ iln(l - 2rcosa + r2) [6{t^ \x\) + e{t - h\) - 6{\x\ - t) - e{\x~ h\ - t)]^ , (9a) 
A\{t,x,0) = -A^(<,x,0) , (9b) 
A\t,x,Q) = 1^(1x1 - t) ln(77>/6) ^ - e(]x-b\ - t) ln(^>/&) 

- i ln(l - 2r cosa + r^) f [0(t - \x\) - 0(|a;| - i)] ^ - [0{t -\x~ h\) - 6{\x - 6| - i)] ^ " ^ 



X X — b 



1 

— arctan 
2 



r cos a 



sin a \ l^l |a; — b\ 



cos a 



(9c) 



The fields Af, A\, and Al {A° = on account of A+ = -A_) are shown for times t/b = 0.3, 1.2, 2.1, and 3 
in Figs. 2-4, respectively. Note that a.tt= \x\ and t = \x_—h\ the fields are continuous, but not continuously 
difFerentiable. This can be most clearly seen in the z-component of the gauge field, cf. Fig. 4. This behavior 
then leads to discontinuities in the electric and magnetic field strengths 

El = -2a_yl^ , (10a) 

E} = -V2d+A^ , (10b) 

Bt = d^Al - dyAl , (10c) 

Bf = V2dyAl , (lOd) 

Bf = -V2d^^Al . (lOe) 

The explicit expressions are lengthy and shall not be given here. Note that the commutator terms in the 
non-Abelian field strengths are absent due to our choice A° ~ S'^^. These terms are in any case of order g"^, 
and thus negligible in the perturbative solution scheme employed here. 
In the field energy density 

7i^i(E? + B?), (11) 

the discontinuities at t = \x\ and t — \x~b\ show up as a series of "peaks" due to the finite resolution of 
the plot routine, cf. Fig. 5. Amusingly, the space-time structure of the field energy density resembles that 
of two expanding "hot spots" in Ref . ! 

Note that to order g^, self-interactions of the produced gluon field are not included. Thus, to this order 
the gluon field produced in an ^ + A-collision is - up to different orientations in color space - a simple 
superposition of the fields generated in collisions of individual charges as shown in Figs. 2-4. Also, to this 
order the field energy density (Fig. 5) is additive. 

The resulting space-time pattern of fields, field strengths, and field energy density can be rather compli- 
cated, depending on the initial positions and color orientations of the colliding valence charges. Also, as was 
pointed out in Ref. |po[| , it is far from homogeneous in the transverse plane in a single event. The subsequent 
time evolution of such inhomogeneous patterns was studied in ||20[] via ideal hydrodynamics, assuming local 
thermalization of matter beginning tq ~ 0.5 fm after impact. This value for the thermalization time scale is 
derived from estimates based on kinetic theory ||2l|,^ . 

Fig. 2: The ^f-component of the radiated gluon field (in units of 1/6) in the transverse plane for times 
t/b = 0.3, 1.2, 2.1, and 3 (a-d). The strength of the color field is A = Sir^/g^. 

Fig. 3: As in Fig. 2, for A\. 

Fig. 4: As in Fig. 2, for Af. 

Fig. 5: As in Fig. 2, for the field energy density. 
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In the remainder of this work we want to address the question of thermalization in the framework of the 
classical approach. Physically, the interactions between the produced gluons, which have been neglected so 
far, will lead to thermalization. These interactions are represented by the non-Abelian terms in the equations 
of motion (H). Since they are of higher order in in the perturbative solution scheme presented above they 
are not considered in the construction of the field to a given order, but only enter as source terms for the field 
to next order in g. Note that "thermalization" refers here to the thermalization of the classical, "soft" gauge 
field degrees of freedom on account of the non-Abelian terms in the field equations, which generate dynamical 
chaos. A priori, that kind of thermalization differs from establishing kinetic equilibrium among "hard" field 
quanta (i.e. particles) through collisions. It also differs from "color" equilibration among particles due to 
color diffusion processes p^ ]. 

In order to address the question of thermalization of the classical gauge fields one would need to self- 
consistently solve the full Yang-Mills equations, including the non-Abelian terms [|^. Because this is a 
numerically rather involved procedure, we here take only a first exploratory step in that direction. 

We want to test the stability of the perturbative solution presented above under (small) perturbations 
obeying the full Yang-Mills equations. If the perturbative solution is stable against these perturbations, the 
latter will either propagate freely through the system, or eventually be damped out. Then, the non-Abelian 
terms are not of major importance for the further time evolution of the system. The perturbative solution, 
which corresponds to non-interacting, free-streaming gluons, dominates the behavior of the system before 
hadronization. 

If, on the other hand, the solution is unstable, the non-Abelian terms are essential for the dynamics of 
the system. Since they represent self-interactions of the field, which eventually lead to thermalization, the 
existence of such instabilities is indicative for the approach to thermal equilibrium. The growth rate of the 
instabilities provides a first order estimate for the thermalization time scale. 

Let us note that the perturbations of the above mentioned stability analysis actually have a physical 
interpretation. Consider the field produced in one single parton-parton collision. Then, the fields produced 
in collisions between other charges can be viewed as perturbations on this field. 



IV. INSTABILITY OF THE CLASSICAL SOLUTION 



A. Linear stability analysis 

To assess whether a given (approximate) solution Aq{x) to the Yang-Mills equations (|^) is stable, we 
decompose the (true) solution A^^{x) as 

A^'ix)^ A<'^{x)+a''{x) , (12) 

where is a small perturbation. With this decomposition, Eqs. (|^) read to linear order in a^: 

□a^ + ig[2 d^Ai^ - ^^'A'(^ ~ ig A^] , a^] - ig [Ao^ , 2 d^a^' ~ d^'a"" ~ ig [A^, a^] + ig [A^, a"] ] 

= - DA^ + ig [Ao, , 2 d'^A^^ - d^Al - ig [A^, Af^] ] = > . (13) 

Since these equations are linear, we may assume without loss of generality that the perturbation is of the 
form 

a''(x) =a(^T°e''=-^ (14) 

(and obtain a general perturbation by linear superposition). With this assumption, the set of partial differ- 
ential equations (^3|) becomes a set of 4 {N^ — 1) algebraic equations for each given wave vector fc^ — (oj, k) 
of the perturbation, 

M^7(fc,Ao)a^ = J,^e-*^ (15) 

where 
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M^,^{k, Ao) = -kHabg'^" + gfabc (2 d'^Aj;^ - 9^^^, - 2ik ■ Aoc g^" + ik^^A^^) 

+g^facefedb [Aoc ' Aod ff'^" " A^^A^J - g^ fabcfcdeA^od^'^^ ' (16) 



Note that M is imphcitly space-time dependent through A^^l^x). For Nc — 2, i.e., fabc = ^abc, Eq. (16) can 
be further simphfied using tabc^cde = SadSbe - ^aehd- 

Kbik, Ao) = -k^abg'"'' + geabc (2 - d'^A'^^ - 2ik ■ A^cg'^'' + zfc^A^J 

+9\Al^Kb - 2Aoa^ofc + ^abA^ocA'-oc + g'-'Aoa ■ A^b - g^^Sab Aoc • Ao.) . (17) 

In order to have non-trivial solutions, we have to require detM(A;,Ao) = 0. This equation has in general 
4(Af^ — 1) complex roots a;,?(k, Ap) which correspond to 8 (iV^ — 1) modes uj.i{\i,Ao) ~ ±yjijjf{\i, Aq) for 
propagation of the perturbation (p^. Note that the dispersion relation Wi(k, Aq) for the perturbation ( p^ ) 
is implicitly space-time dependent through Aq(x). In other words, the perturbation (|lj) is an undistorted 
plane wave only for a constant (i.e., static, homogeneous) field Ag. The space-time symmetries of this 
constant field then refiect in a space-time independent dispersion relation a;i(k, Ao)j and plane waves are 
the true eigenfunctions of the perturbation (p^). 



B. The decay of a constant magnetic field 

As a simple example, let us study the stability of a constant SU(2) magnetic field (see also |^6j), for 
instance (we omit the subscript "0" in the following) 

Bl = H = const. , (18) 

which can be realized by the choice of gauge fields 

Ai = -Hy , aU other A^ ^ . (19) 

In this case, most of the elements of the matrix M are zero, and its determinant is readily evaluated. 
Among the 8{N^ — 1) = 24 roots uiiCkjAo), eight modes simply represent free (harmonic) propagation of 
the perturbation, = |k| , i = 1, ...,4, and uji = — |k| , i = 5, ...,8. In the matrix M, these modes are easily 
identified as corresponding to perturbations which have the same color as the background magnetic field. 
Then, the commutator terms in the equations of motion vanish, rendering them effectively Abelian. The 
Abelian equations of motion are, however, stable under perturbations. 
Furthermore, one finds four stable propagating modes 

c^. = y(fc^^^i/^7TA|Tfc2 , i = 9,10, uo,^-^{k--gHy)-' + kl + kl , i = ll,12, (20) 

and four stable modes 

u^ = y/(fc^ + gHy)^ + kl + kl , 2 = 13, 14 , w,; ^ -y^(fc^ + gHy)"^ + k^ + k'^^ , z = 15, 16 . (21) 

Finally, there are four modes which are solutions of the complex equation 

Lj^-Lu' [2k' -3gHyk- + {gHy)'] 

+ k^ [k^ - 3gHy k^ + (gHy)'] - 2{gH)' [l + ikyy - {k^y)'] - k^{gHy)' = , (22) 

and four modes which obey the same equation with the replacement H —H. Since (^2|) is effectively 
second order in w^, it can be easily solved. 

For our purpose, however, it is sufficient to consider the special case of a non-propagating perturbation, 
k = 0, at position y — 0. Then, both Eq. ( p^ ) and the corresponding one with H —H reduce to 
uj'^ — 2{gH)''. Obviously, the four (two-fold degenerate) roots of these equations are uJi = 2^/^^/gH , i — 
17, 18 , uj., = -2^/^y/gH , i = 19, 20 , and = i 2^1^^^ , i = 21, 22 , tj; = 2^'^^/^ , « = 23, 24. The first 
four modes are just oscillations, while two of the remaining four are exponentially damped perturbations. 
The other two, however, correspond to unstable, exponentiallysrowmg perturbations. Their rate of growth 
is 2^/'^^/gH. We have thus recovered the well-known result |1^ that constant (color-) magnetic fields are 
unstable. Note that the same analysis can be done for constant electric fields, with analogous results. 
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C. Instabilities in the classical gluon field generated by a collision of two color charges 



In this subsection we want to address the question of stabihty of the gluon field radiated in a collision 
between two color charges. Again, we focus on the case of SU(2). Following Subsection IV. A, one would have 
to compute det M{k, Aq) for given fc'^ at any point in space-time. Even for SU(2) this is a formidable task, 
in particular since now, contrary to the case of a constant magnetic field studied in the previous subsection, 
most elements of M are non-zero and, even at z = 0, complicated functions of t, x, cf. Eqs. (pa|-p^) . 

We therefore restrict ourselves to an exemplaric stability analysis at one single point in the transverse 
plane, x = (&/2,0), i.e., half-way between the two colliding charges in the reaction plane. This point is 
obviously a point of special symmetry. Indeed, one finds that all magnetic fields and the y-components of 
the vector potential and the electric field strength vanish here. For times t > b/2, the only non- vanishing 
quantities entering M are: 

^f--0^1n(l + 6V4i^), (23a) 
Al^-^-\n{ ^ , , (23b) 



7r2 t \l + b^/4t^ 

9oAf = ^^-,—^^, (23c) 



doAl 



87r2 i3 1 + ^2/4^2 

87r2 i2 



2 , / 62/^2 

' In ' 



l-t-62/4i2 \l + b^/4:t 



(23d) 



4^N-f^ilnf^^) , (23.) 



7r2 i2 + 62/4^2 

The time evolution of Af, Af, of the electric fields Ef ~ —doAf and = —dzA^ — doAf, and of the field 
energy density H = [{Ef)'^ + {Ef)'^]/2 is shown in Fig. 6. Note that the energy density drops ^ t~^. 

The stability analysis proceeds as in the last subsection. For further simplification we assume a homoge- 
neous, non-propagating perturbation — {uj, 0). The explicit form of the matrix M(fc, Aq) shows that there 
are eight stable modes with = 0, z = 1, 8. At finite k, these would correspond to the freely propagating 
modes found in the presence of a constant magnetic field. Again, the propagation is undisturbed, since the 
perturbations carry the same color as the background field, a — 1. 

We furthermore identify four stable oscillating modes uJi = g^/ (Af )2 -|- (Af )2, i = 9, 10, and uJi = 
—gyJ{AlY + (^i)^i * — 111 12, corresponding to perturbations with colors 2 and 3 and polarization in 
y-direction. The final twelve modes emerge as the complex roots of 



-C^2 



= (24) 

fa2+a2 (ao-iw)a^ 2d^a° + {d^ - iuj)a'' 

-uj^ + al + al -{do-iuj)a' -2d^a^ - {do - iuj)a'- 

29oa"= -uo^ + al -a^a^ 

2800" -uo^ + al -a^a^ 

2900^ + 9^0" -a"=a^ -^2 + a2 q 



-29oa^-92a° -a^a^ -w^ + 

where we have defined = gA'^. Evaluation of the determinant with the help of Mathematica shows that 
the twelve roots are two-fold degenerate solutions of the complex sixth order polynomial 

= (^2 _ a2)2 _,_ ^2 (2 e2 + a,a"9oa^) - (2 a • e -I- a"9^a°) (a • e - a'd,a°) 
+ iuj{uj'^ -a?){2a.-e + a^d^a°) , (25) 

where e = gEi. Decomposing uj — Keuj + iliauj, one ends up with two (real) equations for the real and 
imaginary part of lu. We determined their corresponding roots graphically in the complex cj-plane. Three 
of the six roots have identical imaginary parts, and differ just in the sign of their real part. In Fig. 7 we 
show the time evolution of real and imaginary parts for the three roots with positive real part. 
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FIG. 6. (a) The time evolution of Af (solid) and Ai (dotted) a.t z = 0, x = {h/2, 0). The time is in units of 6, 
the transverse distance between the colliding charges. Fields are given in units of 6"^ (b) The time evolution of the 
corresponding electric fields El (solid) and (dotted), in units of h^^ . (c) The time evolution of the field energy 
density (in units 6~*). For all quantities, the strength of the color field is A = S-ir'^/g^. 
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As one observes in Fig. 7a, one of the imaginary parts is always negative (denoted here by Ima^i), corre- 
sponding to an unstable, exponentially growing perturbation. The other two are either damped (hnuj^ > 0), 
or almost stable (Imw2 — 0). The instability of the first perturbation becomes weaker in time, due to the 
fact that the fields and field strengths decrease (cf. Fig. 6). The rate of growth F = jlmcjil of the instability 
is, however, at least for initial times t ^ b of order b^^. 

The existence of such instabilities and the order of magnitude of their rate of growth is remarkably stable 
against variations in the strength of the color field, A. For instance, at early times t ^ b, the exponential 
growth of the most unstable mode is only about a factor 2 stronger for A being an order of magnitude 
larger, and only about a factor of 3 smaller if A is decrea sed by a n order of magnitude. Since the field 
strengths in the present case are proportional to A, cf. Eqs. (23c-23e), this suggest a power-law dependence 
F ~ A", a ^ 1/2, in rough agreement with the one found in the previous subsection, F ^ ^/H . 



V. CONCLUSIONS 



In this work we have explicitly shown the space-time structure of the classical gluon field produced in a 
collision between two ultrarelativistic color charges. The gluon field radiated in a nuclear collision is, to first 
order, a superposition of such fields. We have demonstrated that this field is unstable against perturbations 
on account of the non-Abclian terms in the Yang-Mills equations. The self-interaction of the produced gluon 
field represented by these terms has been neglected in the perturbative approach of Refs. 10| to construct 
the solution. Our results are in agreement with earlier studies of classical instabilities in non-Abelian gauge 
theories 15 1, and once more emphasize that the transverse field degrees of freedom are of crucial importance 
for the dynamics of high-energy collisions. 

The growth rate of the instability was found to scale with the inverse distance between color charges, 
F ^ b"^ . This implies that, in collisions of large ensembles of color charges (such as nuclear collisions), it 
should scale like the square root of the transverse area density of color charges, F ~ /i, since on the average 
the distance between color charges in the transverse plane decreases inversely proportional to the root of their 
area density This argument is supported by the fact that on the average the gluon field energy density 
(G^) ~ /Lt'* (cf. the calculation of the field energy in fl^), such that F, which was found to be proportional 
to the square root of the field strength, should scale as (G^)^/"* ~ /i. 

We can also find the dependence on the coupling constant. The gluon field strength G is proportional to 
5^, and the results of Section IV. B suggest F ~ \/gG ~ g^^jL. This result is not unexpected, since this is the 
only combination of the coupling constant and the single dimensional scale in the problem, which does not 
contain fi and hence can occur in a classical theory. 

The results of Section IV. C indicate that the coefficient c in the relation F = cg^^ is of order unity, 
such that F ~ (f; ~ 1 for all practical purposes). As has been found in Section IV. C, the uncertainty is 
only a factor of 2, even if the strength of the color field varies by an order of magnitude. Identifying the 
thermalization time scale Tth = F~^, for /x ~ 400 MeV, as estimated by the authors of |^ for RHIC collisions, 
we can expect thermalization times of the order of Tth — 0.5 fm. For LHC, /j, ~ 1 GeV, and Tth — 0.2 fm. This 
is in agreement with estimates of the equilibration time deduced from the maximum Lyapunov exponent in 
random ensembles of classical SU(2) and SU(3) gauge fields Surprisingly, this equilibration time scale 
of classical color fields agrees also with estimates for the time scale for kinetic equilibrium due to collisions 
between "hard" field quanta |^ . Given the fact that the latter results were obtained based on a conceptually 
quite different framework, the agreement is truly remarkable. We conclude that with equilibration time scales 
of this order the odds for producing an equilibrated QGP in RHIC and LHC collisions seem to be good. 
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